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I. INTRODUCTION 

Many modern nanoelectronic devices such as quantum 
dots rely on reduced dimensionality. Two-dimensional 
(2D) electron systems can be well described using a 2D 
Hamiltonian in the effective mass approximation.^* 2 - The 
density-functional theory (DFT) can be applied to de- 
scribe electronic properties in 2D as well as in three di- 
mensions (3D). However, 3D functionals for electronic 
exchange and correlation are not, in general, as such ap- 
plicable for systems with reduced dimensionality^— 

For electronic systems confined in 2D there exists, in 
addition to two local-spin-density approximation (LSDA) 
parametrizations,^ more recent approximations such as 
local functionals a Thomas-Fermi-type explicit func- 
tional of the density^ generalized gradient approxima- 
tions (GGAs)r^ii^ and several Laplacian-level function- 
als (meta-GGAs)r^~— describing either the exchange or 
correlation energy /potential or both. Also the optimized 
effective potential method 2 - and orbital functionals have 
been applied in 2D. 21 However, these approximations and 
their possible combinations remain relatively untested so 
far, and their predictions have not been extensively com- 
pared against one another or against same accurate ref- 
erence data. 

In many cases, approximations for the exchange and 
correlation energy functional are derived assuming that 
the electron density is slowly varying and by start- 
ing from many-body calculations made for homogeneous 
electron gas (for 2D examples, see for example Refs. 
and H) and the local-density or local-spin-density ap- 
proximation (LDA/LSDA). An alternative route taken 
in this work is to consider inhomogeneous model sys- 
tems. Starting from accurate densities and total ener- 
gies, one can obtain accurate exchange and correlation 
potentials and energies for the model systems. Given 
an accurate ground-state density of a model system, the 
Kohn-Sham potential and the exchange-correlation po- 
tential follow uniquely by virtue of the Hohenberg-Kohn 
theorem. Since the exact exchange-correlation functional 
is unknown, the potential can only be obtained by itera- 



tive potential inversion techniques, which can provide the 
potential reproducing a given accurate density. Approxi- 
mate exchange-correlation functional should produce for 
the same density the same exchange and correlation en- 
ergies and potentials so that both the energy and the 
electron density would converge to their correct values. 

The main focus of this Article is on generating accu- 
rate reference data for testing, comparing, and creating 
new functionals and their parametrizations for the elec- 
tronic exchange and correlation in 2D. In order to be able 
to access the exchange-correlation potentials and ener- 
gies of our model systems, we calculate highly converged 
electron densities using the exact diagonalization tech- 
nique. The reference data can be used to benchmark ap- 
proximate exchange-correlation functionals. Most impor- 
tantly, this can include comparing approximate exchange 
and correlation potentials, whose accuracy determines 
the accuracy of the electron density in self-consistent 
DFT calculations, with our accurate ones. To give an 
idea of the accuracy of present functionals, we study the 
accuracy of the LSDA parametrization by Attaccalite et 
al& The LSDA is a unique functional, which becomes ac- 
curate in the limit of slow density variations, and many 
approximations build on top of it. As such, it is the most 
important starting point of comparison. 

The Article is structured as follows. In Sec. HH we 
review briefly the formalism related to DFT calculations 
within the Kohn-Sham method. Section [IIII presents our 
model systems and the specifics of our many-body and 
DFT calculations. Section IIVI contains our results and 
discussion, and Sec. IVl presents our conclusions. 



II. FORMALISM 

In the spin-polarized version of the Kohn-Sham 
method for the density-functional theory, the energy 
functional is written in terms of the two spin densities 
rif(r) and n^(r) and the total density n(r) = n-j-(r) + 
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where the summation in the first term is over occupied 
Kohn-Sham orbitals. Here we use the Hartree atomic 
units. Above, the first term is the kinetic energy of an 
auxiliary non-interacting system with spin densities equal 
to those of the interacting one, the second one the classi- 
cal Hartree energy, the third one the interaction energy of 
the electron system with the external potential, and the 
last term the exchange-correlation energy, which needs 
to be approximated. The orbitals ip? (r) of the auxiliary 
non-interacting system are solved from the single-particle 
equations 



<ff(r) 



C(r)=efC(r), 



(2) 



in which the effective potentials can be written as 

< ff (r) = J dr> + v ext (r) + < c (r), (3) 

where v cxt (r) is the external potential and i£ c (r) the 
exchange-correlation potential, a functional derivative of 
the exchange-correlation energy, 



< c (r) 
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The densities expressed in terms of orbitals are 



K(r)| : 



(4) 
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where one again sums over occupied orbitals. Equa- 
tions (J!])-© are iterated self-consistently until the effec- 
tive potentials v% s (r) and densities n a (r) are consistent 
with one another. 

The Hohenberg-Kohn theorem applied to the non- 
interacting Kohn-Sham system guarantees that if the ex- 
act ground-state densities n a (r) are known, the effective 
potentials v° s (r) follow uniquely, as long as the densities 
are noninteracting v representable. Several numerical al- 
gorithms^— exist for finding the effective potentials cor- 
responding to given ground-state densities. 

In this work, we construct numerically accurate 
exchange-correlation potentials for specific model sys- 
tems in two dimensions according to the above prescrip- 
tion, by starting from accurate wave function calculations 
yielding the electron densities. We compare the obtained 
potentials to ones predicted by the LSDA parametriza- 
tion by Attaccalite et al& and discuss their differences. 



We also obtain information on total energies, exchange 
and correlation energies, and how these are described by 
the LSDA. The data derived from accurate wave function 
calculations can also serve as a database useful for com- 
paring and benchmarking present and future functionals. 



III. COMPUTATIONS AND MODEL SYSTEMS 

In order to be able to analyze the behavior of numer- 
ically accurate exchange-correlation potentials of inho- 
mogeneous electron systems in 2D and compare them 
with ones predicted by the LSDA, we consider a few spe- 
cific model systems and calculate their total energies and 
electron densities using the exact diagonalization tech- 
nique (ED)j2! To reduce the potential inversion proce- 
dure to a one-dimensional problem, we restrict ourselves 
to systems that are radially symmetric. We use either 
a harmonic potential (v ex t(r) = u! 2 r 2 /2) or a circular 
hard- wall potential (infinite potential beyond some given 
radius R and zero within) to confine the electrons. Using 
methods described below, we invert the effective poten- 
tials v° s (r) reproducing the spin densities n a (r) of each 
given model. 

In the special case of two electrons in a spin singlet, 
the effective potential can be obtained directly from the 
density^ The Kohn-Sham orbital is calculated as [see 
Eq. ©] 



V CT (r) = y/n a (r), 



(6) 



and the effective and exchange-correlation potentials can 
then be solved from Eqs. ((2j) and ((3|). 

In a more general case, we use an iterative potential 
inversion algorithm j 25 ' 26 In the course of the iteration, a 
new approximation for the effective potential {v°£ +l (r)] 
is calculated from previous one [ufjj(r)], and the corre- 
sponding density n^.(r) as 



; elf 



r4(r) 



r(r) 



1 <j.if \ 

""off ( r ): 



(7) 



where n a (r) is the reference (ED) density whose generat- 
ing exchange-correlation potential we want to determine, 
and a > is a smoothing parameter removing the effect 
of density tails. In order to keep the iteration stable, we 
start with a large a and decrease its value as the potential 
starts to converge. Also, we take care that the prefactor 
of Eq. (JTJ) does not deviate too much from unity^ The 
behavior of the algorithm has been found to depend on 
the zero level of the potential^ An empirical modifica- 
tion of the scheme and Eq. ([7]) we have found to work well 
for harmonically confined systems is to align v° s (r) to be 
negative and to use an inverse prefactor corresponding to 
iteration to the opposite direction. Then the potential is 
raised where the density is too high and vice versa. 

As we vary the parameters of the confining poten- 
tial and the number of electrons in the system, we 
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get a corresponding set of densities and effective po- 
tentials to analyze. We can, for instance, compare the 
exchange-correlation potentials and energies determined 
by Eqs. © and JTJ to results predicted by the LSDA. 

In the case of the harmonic confinement with varying 
uj, we focus on two electrons in a spin singlet, 4 electrons 
with (L,S) = (0,1), and 6 with (L,S) = (0,0), and in 
the case of the hard- wall confinement with varying R, on 
two electrons in a spin singlet state. For the harmonically 
confined systems we either use expansion in relative co- 
ordinates (two-electron case, number of terms taken high 
enough to give numerically exact results) or the simple- 
harmonic-oscillator basis (4 to 6 electrons) and full ED 
with up to 55 single-particle basis functions. For the 
hard- wall systems, we use Bessel functions and do full 
ED with up to 50 single-particle basis functions. In our 
DFT calculations we use a Bessel function basis. 

The Hamiltonian of interacting Coulomb particles in 
an external potential u ex t(r), 



1 r 



fext(ri) 



-T — 



r; - r, 



(8) 



can be in many interesting examples of confining poten- 
tial expressed via nondimensionalization effectively as 
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where v' cx± is now independent of the parameter charac- 
terizing the confinement whose effect now is to determine 
a unit system, namely natural length and energy scales 
for the specific model potential, and scale the strength 
of the electron-electron interaction. The nondimension- 
alized Eq. © is obtained from Eq. © by substituting 
r — > r/j, choosing 7 2 as the energy unit and the value of 
7 so that r becomes a dimensionless variable, and iden- 
tifying the new v' ext independent of the confinement pa- 
rameter. In Eq. ©, the natural units of energy are j 2 
and those of the length I/7, and the interaction is scaled 
by I/7. For instance, for the harmonic confining po- 
tential Uext(r) = w 2 r 2 /2, 7 = y/uj and v' cxt (r) = r 2 /2, 
and the units for energy and length are uj and 1 / \fu~j 
(HO units). Then in the HO unit system, having a con- 
finement of uj corresponds simply to scaling the interac- 
tion by 1/V^ while keeping the external potential fixed 
(a; = 1). For the hard- wall potential with confinement 
radius R, 7 = 1/R and the units for energy and length 
are 1/i? 2 and R (HW units). 



IV. RESULTS AND DISCUSSION 
A. Energies 

Our model systems are summarized in Table U where 
we also list our ED total energies and DFT total ener- 
gies calculated using the LSDA parametrization by At- 



taccalite et al& All the LSDA results in this article corre- 
spond to ground-state densities of self-consistent LSDA 
calculations. Using ED densities and corresponding or- 
bitals would not affect our conclusions. The potential 
inversion procedure provides us, in addition to the ac- 
curate effective potential, also with accurate Kohn-Sham 
orbitals. Using these and the Kohn-Sham energy func- 
tional [Eq. ([T])] we can calculate the exact ground-state 
exchange-correlation energy, E xc , for any given model 
system. These, along with the corresponding approxi- 
mate LSDA values, are also listed in Table |U In the case 
of the spin-singlet two-electron systems the exchange en- 
ergy is simply 28 



E x [n t , ni ] = -- 



dr dr' 



, n(r)n(r') 



(10) 



i.e., minus one half times the Hartree energy. This allows 
us to easily decompose the exchange-correlation energies 
of these systems into exchange and correlation parts. 

Some ED reference results exist for the total energies of 
harmonically confined systems in the literature. For the 
two-electron systems our energies are lower than results 
of Ref. 29:, which are 3.72143 (uj = 0.25) and 3.00097 
(a; = 1). For uj = 1 there exists an analytic solution 
with an energy of 3 (Ref. [30h . For the 4-electron system 
with uj = 0.25 our energy is only slightly higher than one 
calculated by Mikhailov^i who obtained 13.6180 using a 
larger one-body basis. For 6 electrons with (L, S) — (0, 0) 
and uj = 0.25 Rontani et ai— have obtained the energy 
of 27.98 using 36 single-particle states, which is higher 
than our result calculated with 55 states. 

The ED and LSDA total energies are rather consis- 
tent for all systems, the latter one being curiously consis- 
tently higher, despite the fact that the DFT total energy 
is not guaranteed to be variational once the exchange- 
correlation energy functional is approximated. Similarly, 
the LSDA exchange-correlation energies are higher, the 
exchange component being clearly underestimated. The 
lower correlation energy predicted by the LSDA partly 
compensates this. This finding of error cancellation is 
consistent with the results for He isoelectronic series in 
3Di2£ The mechanisms behind the underestimation of the 
exchange energy and cancellation of errors between ex- 
change and correlation energies are the same for the 2D 
as for the 3D LSDA. The underestimation of exchange 
energy is due to self-interaction, whereas the compen- 
sating effect of the overestimation of the correlation en- 
ergy arises from the exchange-correlation hole sum rule. 
Since the 2D LSDA is based on a physical system, the 
2D uniform electron gas, the sum rule is fulfilled and the 
integrated errors of the exchange and correlation holes 
canceL 33 ' 34 

Figure [1] shows a graphical representation of LSDA's 
relative errors in the total energy and the exchange- 
correlation, exchange and correlation energies for the 
two-electron systems in harmonic and hard-wall confine- 
ments as a function of the characteristic length scale of 
the system, 1/y/uJ or R. The relative errors for the total 
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TABLE I. Summary of our model systems and calculated energies. The table includes the systems' descriptions, shape of 
the confining potential and the related parameter R or our ED total energies, i?toti the LSDA total energies £^t^^, 
calculated using the parametrization of Ref. [s| the exchange-correlation energies inverted from the ED calculation, _E XC , and 
the corresponding LSDA values, £jf DA . For two-electron systems we also show their decompositions into exchange and 
correlation energies, _E X and E c . Natural units determined by the external potential are used throughout. 
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and exchange energies are rather constant as a function 
of the systems' length scale. The same applies to the 
exchange-correlation energy, which consists mostly of the 
exchange one. The relative error in the correlation energy 
is larger at stronger confinements (large ui j small R) and 
becomes smaller at weaker confinement. This does not af- 
fect the overall picture much since the correlation energy 
is a small fraction of the exchange-correlation one, espe- 
cially at the weakly interacting (large u> / small R) limit. 
The strongly confined systems are less uniform than their 
counterparts in weaker confinement. Therefore, the ac- 
curacy of the LSDA is worse in this limit. The accuracy 
of the LSDA energies is improved with increasing par- 
ticle number. For two electrons, the level of accuracy is 
the same for both the hard-wall and harmonic potentials. 
For 4 electrons, the relative error in the total energy is at 
most 0.6% and in the exchange-correlation energy 2.4%, 
and for 6 electrons 0.3% and 1.3%, respectively. 



B. Densities and potentials 

We begin the comparison of exchange-correlation po- 
tentials from the two-clcctron systems in spin singlet. 
Figure [2] shows electron densities and different potential 
terms (exchange and correlation, exchange, correlation) 
for the harmonically confined electrons with varying u> 
and hard-wall systems with varying R. The scaling we 
use when representing the data is described and moti- 
vated below. For the two-electron systems in spin-singlet, 
the exchange potential is2£ 



0.2 | | i | , | i | , | r 



0.1 - ♦ — — ♦ 




123456789 10 
R (a.u.) 



FIG. 1. (Color online) Relative errors in the LSDA energy 
terms (total, exchange-correlation, exchange, and correlation 
energies) for (a) two-electron systems in harmonic confine- 
ment and (b) two electrons in a hard-wall confinement as a 
function of the characteristic length scale of the system, 1/%/^ 
or R. 



When visualizing the potential terms we use natural 
units and scale the exchange potential by 7 and the cor- 
relation one by 7 2 . Since the exchange-correlation po- 
tential consists mostly of the exchange one, we scale it 
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FIG. 2. (Color online) Densities (a), exchange-correlation potentials (b), exchange potentials (c), and correlation potentials (d) 
as a function of the distance r for two electrons in a spin singlet in harmonic confinement u, and (e)-(h) the same quantities 
for two-electrons in the hard-wall confinement with radius R. For the latter the results for R = 2, . . . , 10 are shown. Solid lines 
are the accurate results and dash lines the LSDA ones. Thin dotted lines in (b), (c), (f), and (g) show the exact asymptotic 
—1/r behavior of the exchange potential. Natural units determined by the external potential are used throughout. 



similarly. These scalings provide energy scales at which 
the potentials are of comparable magnitude and their fea- 
tures can be easily compared. This can be explained 
using the following scaling relations for the exchange en- 
ergy^ 



and the correlation energy^ 

E c [n x ] = A 2 £ C V >]- (13) 

Above, A is an arbitrary scaling parameter not necessarily 
referring to a transformation between two unit systems, 
and the scaled density ri\ is denned as 



E x [n\] = XE x [n], 



(12) 



n A (r) = A d n(Ar), 



(14) 
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where d is the dimension (here 2), and Ej [n] is the 
density functional for the correlation energy for a sys- 
tem with density n but with electron-electron interaction 
scaled by 1/A. The corresponding scaling relations for the 
exchange and correlation potentials are analogous. For 
the exchange potential- 
ly (fax]; r) = \v x ([n]; Ar), (15) 
and for the correlation potential, 

v c ([n x }:r) = X 2 v 1 J x ([n];Xv). (16) 

The latter one follows from Eq. (|TU1) similarly as Eq. (fT"5)) 
is derived from Eq. ([H]) (see Ref.]33). We apply Eqs. (15]) 
and (fl"6|) in such a way that n corresponds to the density 
of a given system expressed in atomic units and n\ the 
same density scaled to natural units (A = 7). Obviously 
the interaction strengths of the functionals in Eqs. (fT5")) 
and (|16p do not match with those of the above unit sys- 
tems, 1 and I/7, [see Eqs. © and ©], but, nevertheless, 
the scaling relations motivate a consistent visual repre- 
sentation. 

There is quite a good agreement in the densities, 
Fig. Ufa), between the accurate ED calculations and the 
LSDA ones for the harmonically confined two-electron 
systems in the weakly correlated cases (large w). At 
smaller uj, the LSDA densities are monotonous while the 
accurate ones develop a side peak. Also for the electrons 
in the hard- wall trap [Fig. (He)], the agreement is best in 
the weakly correlated limit (small R). The trend in the 
accuracy is here opposite to the one seen for exchange 
and correlation energies in Fig. [TJ where the energies are 
in better agreement in the strongly correlated (uniform 
system) limits. Integrated quantities such as the energy 
can behave differently from local quantities such as the 
density or the potential in the sense that integrated quan- 
tities can be more accurate due to cancellation of local 
errors. Especially in the case of the harmonic potential 
the agreement in the densities is bad only at a small area 
at the center of the quantum well. 

In the exchange-correlation potential of the harmoni- 
cally confined electrons, Fig. [U[b), the most marked dif- 
ferences are the different asymptotic behavior (— 1/r vs 
Gaussian decay) due to the lack of self-interaction cor- 
rection in the LSDA exchange potential, see also the 
exchange potential in Fig. dfc), and the resulting ver- 
tical shift. For the hard- wall systems, the differing be- 
havior of the potentials close to the wall is even more 
pronounced due to the rapid decay of the charge den- 
sity, see Figs. HJf)-(g). The shapes of the exchange- 
correlation potentials close to the potential well center 
differ, especially in the case of harmonically confined 
systems. The unphysical monotonousness of the LSDA 
exchange-correlation potentials of the harmonically con- 
fined systems arises from the LSDA's local character and 
the monotonic behavior of the LSDA densities. Separate 
comparisons of the exchange and correlation potentials of 
these systems, Figs. [UJc) and^d), reveal that the non- 



monotonicity of the accurate exchange-correlation poten- 
tial arises from the correlation part. It is also noteworthy 
that the correlation potential changes sign to positive at 
large radii and decays to zero from above. This behav- 
ior is familiar from 3D systems! 23 ' 28 In the hard- wall sys- 
tems, [Figs.[Uf)-(h)] the accurate correlation potential is 
again non-monotonous making the exchange-correlation 
potential non-monotonous for all values of R. On the 
other hand, the behavior of the accurate exchange poten- 
tial is, for almost all values of R apart from the largest 
considered, monotonous, whereas the LSDA exchange 
potential displays a stronger non-monotonous behavior 
following the trends of the local electron density. In con- 
clusion, for the hard- wall system the LSDA exchange and 
correlation potentials take different roles than the ac- 
curate ones and the resulting canceling of errors leads 
to exchange-correlation potentials that are qualitatively 
correct. This is similar to the error cancellation between 
exchange and correlation energies discussed above. 

The next closed-shell system is 4 electrons with 
{L,S) — (0,1). Figure [3] shows the spin densities and 
spin-dependent exchange-correlation potentials for vary- 
ing u). For the majority spin [Fig.[3][a)-(b)], both the den- 
sities and exchange-correlation potentials are described 
rather accurately. For the minority spin [Fig. EJc)-(d)], 
the accuracy decreases with decreasing ui. The LSDA 
exchange-correlation potentials are, however, qualita- 
tively correct apart their wrong asymptotics and less re- 
pulsive shape at the origin. 

Finally, Fig. 0Ja) shows densities and Fig. Sfb) 
exchange-correlation potentials for 6 electrons with 
(L, S) = (0, 0) in harmonic potential with varying uj. For 
these unpolarized systems with higher average electron 
densities, the accuracy of the LSDA is remarkable both in 
the densities and exchange-correlation potentials. Only 
at small lo the LSDA and accurate exchange-correlation 
potentials start to display differences. 

In addition to plotting the model systems' exchange 
and correlation potentials as a function of the position we 
can analyze how they look when shown as a function of 
the local electron density and compare them in this repre- 
sentation with the LSDA. Deviations from the LSDA give 
a measure of the non-locality of the exchange and corre- 
lation effects. Figure [5] shows a comparison for our un- 
polarized systems with 2 or 6 electrons. Figures EJa)-(b) 
show the exchange-correlation and correlation potentials 
of the harmonically confined two-electron systems, both 
accurate and LSDA ones, of Fig. (5Jb) and (d) as a func- 
tion of the local electron density, and Figs.[5jb)-(d) those 
of the hard- wall systems of Figs. [5Jf) and (h). Figure^e) 
shows the exchange-correlation potentials (see Fig. [4j of 
the 6-electron (L,S) = (0,0) systems as a function of 
the local density. When looking at the total exchange- 
correlation potential, the accurate result and the LSDA 
one behave, in general, similarly. The vertical shift due to 
differing asymptotic behavior of the exchange potential 
is apparent. In the cases in which the accurate density 
is non-monotonous it becomes clear that the exchange- 




FIG. 3. (Color online) The majority spin density (a) and majority spin exchange-correlation potential (b) as a function of the 
radius for 4 electrons with (L,S) = (0, 1) in harmonic confinement w, and (c)-(d) the same quantities for the minority spin. 
Solid lines are the accurate results and dash lines the LSDA ones. The thin dotted line shows the exact asymptotic —1/r 
behavior of the exchange potential. HO units are used. 



correlation potential cannot be expressed as a function 
of the local electron density but it is a non-local func- 
tional of the density. In the case of the 2-electron sys- 
tems, we can identify the correlation potential, Figs.tSfb) 
and EJd) as the main-source of the non-locality close to 
the density maxima. Other manifestations of the non- 
locality include the above-mentioned —1/r asymptotics 
of the exchange potential absent from the LSDA one, 
see Figs. [U[c) and (g), and the differing character of the 
exchange potential of the hard- wall systems, Fig. [Hg), 
see the discussion above. The behavior of the exchange- 
correlation potential as a function of the local density is 
similar between the 2 and 6-electron systems. Also in this 
comparison, the LSDA seems to work better for larger 
particle numbers. If Fig. [5] were plotted using atomic 
units, the LSDA potentials would collapse on top of each 
other. The universality seen in the accurate results of 
Figs OHa), [H[c), and[SJe) is an essentially non-local effect. 
It arises from the correct — 1/r asymptotics of the ex- 
change potentials at the large r (low local density) regime 
of the model systems. 



V. CONCLUSIONS 

In this work, we have extracted accurate exchange- 
correlation energies and potentials for inhomogeneous 
model electron systems in two dimensions starting from 
total energies and electron densities calculated by exact 
diagonalization. We have considered two electrons in a 
spin singlet in harmonic and circular hard-wall confine- 
ments and 4 electrons with (L, S) — (0, 1) and 6 electrons 
with (L,S) = (0,0) in harmonic confinement. We have 
compared our results against results calculated by the 
local-spin-density approximation (LSDA) parametriza- 
tion by Attaccalite et alM- The LSDA appears to describe 
these systems relatively accurately. 

Total energies of our model systems are curiously con- 
sistently overestimated by the LSDA. The exchange en- 
ergy predicted by the LSDA is too high but this is par- 
tially compensated by the too low LSDA correlation en- 
ergy. Considering the exchange and correlation poten- 
tials determining the accuracy of electron densities in 
density-functional calculations, the LSDA exchange po- 



8 



0.5 




r (HO units) 



FIG. 4. (Color online) The electron density (a) and potential 
(b) as a function of the radius for 6 electrons with (L, S) = 
(0,0) in harmonic confinement cj. Solid lines are the accurate 
results and dash lines the LSDA ones. The thin dotted line 
shows the exact asymptotic —1/r behavior of the exchange 
potential. HO units are used. 



tential is obviously lacking the exact —1/r asymptotic 
tail because its lack of the self-interaction correction, a 
consequence of the locality of the LSDA. Nevertheless, 
the shape of the exchange-correlation potential is, apart 
from the two-electron systems in harmonic confinement, 
qualitatively correct, sometimes owing to error cancella- 
tion between the exchange and correlation components. 
In general, the LSDA as parametrized by Attaccalite et 
al£- is quite accurate but it is clearly useful to go beyond 
the LSDA and further develop and test semi-local and 
non-local density and orbital functionals for electronic 
exchange and correlation in two-dimensions. Our results 
shown in this Article can act as benchmark data in creat- 
ing, testing and parametrizing exchange and correlation 
functionals for two-dimensional electronic systems. 
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6 electrons with (L,S) = (0,0) in harmonic confinement (see Fig. Solid lines are the accurate results and dash lines the 
LSDA ones. Natural units determined by the external potential are used throughout. 
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